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Abstract 



Based on an exact trace formula for a one-dimensional ray-splitting system, 

we derive novel combinatorial identities for cyclic binary sequences (Polya 

CN . necklaces). 
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I. INTRODUCTION 

Wave propagation in systems with sharp interfaces is a fundamental problem in the nat- 
ural sciences and engineering. Well-known examples include light waves impinging on a 
water-air interface or sound waves propagating in layered media [0] . All these systems have 
one feature in common: the splitting of the incident wave into reflected and transmitted 
components. In the geometrical optics limit of small wave lengths, incident, reflected and 
transmitted waves are described by rays. The rays are split at the interface; hence the name 
"ray-splitting systems" H for the whole class of wave systems with sharp interfaces. Ray- 
splitting systems have recently attracted attention in the context of acoustic, quantum and 
electromagnetic wave chaos [PHT0||. It was shown that the mere presence of a ray-splitting 
boundary can drive an otherwise regular system into chaos 0-ffl. It was also shown that 



ray-splitting systems produce corrections to the Weyl formula |TT| , p!2[| for the average den- 
sity of states that can be computed analytically [0,0. The most conspicuous consequence 
of ray splitting is the existence of non-Newtonian periodic orbits in a ray-splitting system 
that contribute substantially to the fluctuating part of the level density @||.|7]^9|] . The 
existence of non-Newtonian orbits in a dielectric-loaded Bunimovich ray-splitting stadium 
was demonstrated experimentally 0HH • in addition it has been shown recently that ex- 
act trace formulae exist for a class of one-dimensional ray-splitting systems PJTt]|]. These 
formulae were derived using results from quantum graph theory |13|,14|. Considering the 



two-point correlation function of the spectra of special quantum graphs, Schanz and Smi- 
lansky were able to derive novel combinatorial identities [ 15| . This was possible by deriving 



the two-point correlation function in two independent ways, (i) directly using input from 
the quantum spectrum and (ii) using the exactness of the trace formula. Motivated by 
the methods of Kottos, Schanz and Smilansky |T3|-|T5| we show that novel combinatorial 
identities for binary Polya necklaces |16|JT7[] are obtained directly by comparing the spectral 
density of analytically solvable quantum graphs with their exact periodic orbit expansion. 
The plan of this paper is as follows: in section II we present our model system, a one- 



dimensional ray-splitting system, whose spectrum can be obtained analytically. We present 
an exact periodic-orbit expansion of its level density expressed as a generalized Fourier 
sum over binary Polya necklaces. In section III we outline our method for obtaining exact 
combinatorial identities for binary necklaces derived from analytically solvable cases of the 
one-dimensional ray-splitting system. In section IV we present two worked examples that 
yield two infinite sets of combinatorial identities. In section V we discuss our results and 
conclude the paper. 

II. SPECTRUM 

Denote by E n the spectrum of the one-dimensional scaled Schrodinger equation 

-iP"(x) + V x (E,x)ip{x) = E^{x), (1) 

where 

'0, for < x < a, 

V X (E, x) = | XE, for a < x < 1, (2) 

k oo, for x ^ (0, 1) 
and ip(x) = for x ^ (0,1). Scaling potentials of the form (|j) arise naturally in many 
ray-splitting systems, for instance in dielectric- loaded cavities H^J^^H- ^ n this P a P er we 
concentrate on the case E > V\(E,x) for all x. Then, without restriction of generality, the 
scaling constant A can be assumed to satisfy < A < 1. Define E = k 2 , then the spectrum 
of (D) is determined by 

sin [k(a L + a R )} - r sin [k(a L - a R )] = 0, (3) 



where o^ = a,a R = 77(1 — 0), rj = \/l — A, and r = (1— i])/(l+i]) is the reflection coefficient. 



For the derivations below it is useful to define the transmission coefficient t = y/1 — r 2 . 
Therefore, the reflection and transmission coefficients satisfy the relation 

r 2 + t 2 = l. (4) 



Since in this paper we focus on the case < A < 1, both r and t are real and positive 
and range between and 1. Possible quantum phases incurred at reflection or transmission 
events are treated explicitly and separately (see, e.g., (|) below). They are not included in 
r or t. 

We now discuss an alternative method of solving the quantum dynamics in the potential 
(@). This method is based on coding the periodic orbits with the help of symbol strings. We 
denote a bounce off x = by the letter £ and a bounce off x = 1 by the letter 1Z. Words 
formed with these two letters code for periodic orbits in (§). The word £, for instance, 
codes for the non-Newtonian orbit that bounces between x = and x = a (above-barrier 
reflection orbit). The word £TZ codes for the (Newtonian) orbit that bounces between x = 
and x = 1. Since a periodic orbit represented by the word w cycles through the letters of w 
without a well-defined beginning or end, two words w and w' are equivalent in our context, 
and code for the same periodic orbit, if they are of the same length (i.e. they consist of the 
same number of symbols) and their respective symbol sequences are identical up to cyclic 
permutations. Sequences of objects that are identical up to cyclic permutations are called 
(Polya) necklaces p^JT7| . If the number of objects they consist of is two, they are called 
binary necklaces. Apparently therefore, the periodic orbits of (0) can be coded with the help 
of binary necklaces over the symbols £ and TZ. It is remarkable that for @ every Newtonian 
or non- Newtonian periodic orbit can be mapped one-to-one onto a binary necklace. In other 
words, "pruning" is not necessary for the binary necklaces relevant for (0). Every binary 
necklace defines a possible periodic orbit of (0) and vice versa. 

Given two letters, for instance £ and TZ, we can form 2 e words of length £. But, in 
general, many of these words will be cyclically equivalent, and correspond to the same 
necklace. So, how many necklaces of length £ are there? This question is answered by the 
following formula. There are exactly p 



N(£) = ~ £ *(») 2 £/n (5) 

1 n\£ 

binary necklaces of length £, where the symbol Li n\£" denotes "n is a divisor of £" , and 4>(n) 



is Euler's totient function denned as the number of positive integers smaller than n and 
relatively prime to n with 0(1) = 1 as a useful convention. Thus the first four totients are 
given by 0(1) = 1, 0(2) = 1, 0(3) = 2 and 0(4) = 2. We illustrate ([5]) with two examples for 
£ = 1 and 1=1. There are two necklaces for £ = 1, C and TZ. Applying ([5]) to this problem, 
we verify N(l) = 0(1) x 2 = 2. There are three necklaces of length 2, ££, CTZ and 7ZTZ; 
again verified by (|), N(2) = [0(1) x 4 + 0(2) x 2]/2 = 3. 

Given a binary necklace w, we define the following integer- valued functions on w: n-ji{w) 
counts the number of IZs in w, nciw) counts the number of £s, n{w) = nciw) + n-ji(w), 
x{w) is the sum of n{w) and the number of 7?.-pairs in w, a(w) counts all occurrences of 
7£-pairs or £-pairs, (3{w) counts all occurrences of 1ZC or CR and 7(u>) is defined as j(w) = 
2n c {w) + n-Tiiw). Note that the counting of 7£-pairs, £-pairs, CTZ- or 7££-combinations is 
to be understood cyclically, i.e., for example, a (1Z) = 1 and (3{C1Z) = 2. Next we define the 
set W p of prime necklaces as the ones that cannot be written as a periodic concatenation 
of substrings. As shown recently ||10|| , there exists an exact periodic orbit expansion for the 



spectral density of (|I|) in terms of prime binary necklaces, 

•I OO || 

p{k)=p+±- y S w Y [(-l)xM r «W^Ml M e iuS - k , k>0, (6) 

/IT L 

z,/l wGW p >>=-<*> 

where 

S w = 2[n n (w)(j R + n c (w)a L ] (7) 

is the action of the primitive periodic orbit coded by the prime binary necklace w and 
p = (cl + ctjiJ/tt is the average level density. 

III. METHOD 

For special values of the parameters of the potential well @ it is possible to solve 
(|) analytically, thus obtaining directly the density of states p(k). Equating the explicit 
expression for p(k) with the necklace expansion (^), one obtains combinatorial identities for 



binary necklaces. An illustrative example is the case a = 1, for which the spectral density 
of fll]) is given by 

oo 

p{k) = Y, S{k-nm). (8) 



m=— oo 



In this case there exists only one primitive necklace, CTZ, and the necklace expansion of (^) 
yields 

1 oo 

P(k) = - E ***■ (9) 

v=— oo 

Equating (|8]) and (|9|) yields the well-known Poisson formula. 

A comment is in order here. For every finite a the necklace expansion (|6]) involves an 
infinite sum over prime periodic necklaces. At a = 1 this sum collapses to a single term. 
One may ask the question how this singular limit arises. The answer is the following. For 
every finite a the actions of the right-hand lobes of the periodic orbits of (|P is finite. At 
a = 1, these actions are zero. Consequently, all necklaces that represent different prime 
periodic orbits for a^ 1 become repetitions of the Newtonian periodic orbit £TZ at a = 1. 
This is the reason for the existence of only a single prime periodic orbit {CTZ) at a = 1. 

Apart from trivial and well-known identities such as (|8|) and (|9p above, (|6p is a rich 
source of new and nontrivial combinatorial identities for binary necklaces. Specific examples 
are discussed in Section III. Here we outline the general method. 

Equation (H) can be written as 

sinfwifc) — rsm(u2k) = 0, (10) 

where LO\ = cr L + cr R and cu 2 = o L — ctr. A negative u 2 corresponds to a mirror reflection 
of V\(E,x) with respect to x — 1/2. Thus, because of ol^r > 0, and without loss of 
generality, we may assume oj\ > uo 2 > 0. 

In case uj\ and u>2 are rationally related, i.e. uj\/u)2 = p/q, p > q E N and p, q relatively 
prime, ( |10|) is reduced to the algebraic equation 

sin(pujk) — rsm(quk) = 0, (11) 



where u>i = puo and U2 = qui. Using the formula 

sin(rax) = sin(x) f/ n _i(cos(x)), (12) 

where U n -i(x) is the Chebyshev polynomial of the second kind, one obtains 

sin(o;fc) [U p _i(cosu)k) — rU q ^i(cosujk)] = 0. (13) 



It follows immediately from (|13D that in the case of rationally related iO\ and U2 there always 
exists a sequence of roots k^ = irn/u. The remaining roots are determined by 

U^x) - rU q -!{x) = 0, (14) 



where x = cos(uk). Since every root Xj of (|14D gives rise to a periodic sequence of eigenvalues, 



cos(uk^) = Xj, j = 1,2, ...,p — 1, together with the sequence k!£> the spectrum of (pT0|) 

consists of p (possibly degenerate) periodic sequences of roots. Whenever (|14"D can be solved 

analytically, the density of states 

p— 1 00 
P(k) = T, E 6{k-W) (15) 

is known explicitly and together with (^j leads to a host of combinatorial identities for binary 
necklaces. Two examples are presented in the following section. 

IV. COMBINATORIAL IDENTITIES 
Example 1. For ol = <jr equation (|3]) becomes 

sm{2ka) = 0, (16) 



with the solutions k n = nn/ (2a). Note that there is no r-dependence in (|I6"D. The density 
of states is given by 



According to (|7|) and due to ai = <Jr, S w depends only on the binary length of w and is 
given by S w = 2an{w). Thus the sum (||) can be written as 



2a a °° 

p( k ) = i + - E E 

7T 7T 



n(w} 



' _\\xiw) r a(w)^(w) 



Jlimka 



(18) 



in | 



where W n denotes the set of all length-n binary necklaces w, w is the shortest primitive 



code-piece in w and v is the number of its repetitions in w. Comparing the series (|17D and 
(|18|), we see that odd-length and even-length binary necklaces satisfy the sum rules 



J2 n{w) [(-l)x(«VM^(«0l " = 0, m = 1, 2, 

W-v£W2m-l 



(19) 



and 



E n ( w ) 

W-V&V2m 



f_l\x( w ) r a ( w )f^(.w) 



1, m = l,2, 



(20) 



At first glance it may seem surprising that (|20|) is a constant for all r. The solution lies in 
the relation (|]). When properly ordered according to powers of r and t, it turns out that 
(i) individual terms in (|20| ) are of the form r 2p t 2q , where p + q = m and (ii) the coefficients 
in front of the term r 2p t 2q in (^) turn out to be binomial coefficients. Thus, for given m, 
the left-hand side of (|20| ) reduces to (r 2 + t 2 ) m , which, according to (Q) is equal to 1 for any 
choice of r. This explains why the seemingly variable left-hand side of ( PD| ) is nevertheless a 



constant. Thus we obtain from (|20D the following infinite set of combinatorial identities for 
even-length binary necklaces 



1 



- £ nH(-ir*^VaWA. 



m 
s 



s — 0, ..., m, m = 1, 2, 



(21) 



where 5jj is the Kronecker symbol. 

In order to illustrate ( |21| ) let us first focus on the case m — 1. According to (|5|) there 
are exactly three cyclically non-equivalent necklaces of binary length 2 given by w x = CC, 
W2 = CTZ, W3 = 7^.7?.. The necklaces wi and W3 are not primitive. The necklace wi is a two- 
fold repetition of the primitive necklace Wi = C. Thus v\ = 2. An analogous consideration 
for W3 yields w 3 = TZ and z/3 = 2. The necklace W2 is primitive. Therefore w 2 = W2 = CTZ 
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and v-i = 1. The three necklaces Wj, j = 1,2,3, are listed in Table I. Also listed are their 
primitives Wj, the repetition indices Uj, and the values of the functions n(wj), a(wj), /3(wj), 
j(wj) and xiwj). 

We are now ready to check (pT|). For m = 1 we have two choices for s: s = and s — 1. 
For s = we have to scan the three words Wj, j = 1, 2, 3, for Uja(wj)/2 = s = 0. According 
to the entries in Table I, only W2 qualifies and the sum on the left-hand side of (^l|) reduces 
to the single term 

l -n{w 2 ){-iy^ w ^ = 1=( J . (22) 

This shows that ( pT|) is indeed true for the simplest case m = 1, s = 0. For the case m — 1, 
s = 1 we have to check Table I for occurrences with Vj(x{wj)/2 = 1. This is fulfilled for the 
necklaces w 1 and w 3 . We obtain 

' r n( Wl )(-l) Ul ^ Wl) + n(w 3 )(-l) U3x( - W3) ] =1=( j . (23) 



This shows that (|21|) also works for m = 1, s = 1. 

Testing (^l]) for m = 2 involves finding all non-equivalent necklaces of binary length 4. 
According to (|5|) there are exactly six. All six necklaces are listed in Table II together with 
their properties. For m = 2 we have three possibilities for s: s = 0, 1, 2. For s = we have 

to check Table II for necklaces that fulfill i>ja(wj)/2 = 0. Only W4 qualifies. We obtain 

'2' 



n(u>4)/2 = 1, which equals ( | , the binomial coefficient on the right-hand side of (0). For 

,0, 
s = 1 we have to check for Uja(wj)/2 = 1. We find three candidates: w 2 , w 3 and w 5 . This 

time we have to be careful when summing the three terms on the left-hand side of ([H]), since 

z/ 3 x(w 3 ) = 5. Therefore the second term in (^1|) contributes with a minus sign. We obtain 

[n(w 2 ) — n(w 3 ) + n(w 5 )]/2 = 2, which equals , the corresponding binomial coefficient 

w 

on the right-hand side of (^1|). Two necklaces, wi and W6, contribute in the case s = 2 and 

again satisfy ([21]) . 

Example 2. Suppose now that a^ = 2a R . In this case (^ becomes 

sin(A;a/2) Ucos 2 (A;a/2) - r - ll =0. (24) 



This equation has three sets of solutions 



fc<?> = ^arccos(p) + , j = -1,0,1, (25) 

Uj (Jj 



where if = \/l + r/2. The density of states is 



7=-i m=-oo V a a / 



1 oo oo 

— E E e in[afe+2j ' arccos(v)1 = — £ e inafc [2T 2n (<p) + 1] , (26) 

^ 7r j=-l n=-oo ^7T n=-oo 



where T n (x) = cos(narccosx) are the Chebyshev polynomials of the first kind |18|| . Equating 
(EH) order-by-order with the necklace expansion (pi) we obtain the sum rules 



_l\xM r a(w)^ _ r 2\f3(w)/2 



"vy(w),m 



E EtH 

m 2m(-\V /2m - j 
1 + E^ZT Plfrr- m=1.2 127) 



We used formula 22:6:1 of ref. [jn| for the Chebyshev polynomials in (p6[). 



Ordering (p7|) according to powers of r, ([27]) can be reformulated as a combinatorial 
theorem on the set of binary necklaces, in which L beads weigh twice as much as 1Z beads: 

27(10) 

2^7) 

™^2m(-lV' /2m- A (m-j\ 

The condition C in the sum (p8|) is C = 7(u>)|m As — ma(u>)/7(-u;) even. The sum on the 
left-hand side of (p8|) may be empty. In this case the sum is defined to be zero. 

Let us check (|28| ) with the help of a few examples. First we focus on the case m = 1, 
s = 0. In order to fulfill the first part of the condition C in (|28|) we need 7 = 1. This, in turn, 
requires to find a necklace with ri£ = and n^ = 1. There is just one such necklace, namely 
1Z. But it does not fulfill the second part of C. Therefore, the sum on the left-hand side of 
(p8|) is empty, and the left-hand side is zero. The right-hand side adds up to 1 + 1 — 2 = 
and confirms (^) for this special case. For s = 1 we find again that 1Z is the only choice for 
w. But this time the second part of C is fulfilled and the left-hand side of (|28|) is 

10 



7 (7^V_l)Px(W)47(«)-a(W)]/[27(W)] [ V L =1 ( 2 9) 

V[7(^)-a(^)]/[2 7 (^)]/ 
We used a(K) = 1, /?(?£) = 0, 7(7^) = 1 and x(X) =2. For m = 1, s = 1 the right-hand 
side of ( pffi ) consists of just one term, which turns out to be 1 as well. Thus we checked that 
(p8|) works for m = 1. With the help of Tables I and II, other special cases may be checked 
as well. 



V. DISCUSSION AND CONCLUSIONS 

The work presented here is closely related to the theory of quantum graphs |13|-[T5|. 
While the quantum graphs considered by Kottos, Schanz and Smilansky [p^ -[T5| correspond 
to the case of zero potential on the bonds and delta potentials on the vertices, the step 
potentials considered in this paper correspond to constant potentials on the bonds and 
potential steps at the vertices. Thus, although the methods employed in this paper are 
essentially those used previously by Kottos, Schanz and Smilansky, we obtain a different 
class of combinatorial identities that apply to cyclic binary codes (Polya necklaces). Another 
difference concerns the derivation of identities. While Kottos, Schanz and Smilansky use 
a route that involves two-point correlation functions, we show that novel combinatorial 
identities can be obtained directly from the periodic orbit expansions of explicitly solvable 
cases. These minor differences notwithstanding the central idea for generating entirely new 
classes of combinatorial identities is the same: Combinatorial identities can be obtained 
whenever a quantum system admits of (i) an explicit analytical solution and (ii) an exact 
periodic orbit expansion. 

In addition to the two examples presented above, there exist many other cases in which 
(|Tjj) can be reduced to a low-order polynomial that can be solved by elementary means. 
Examples are the cases p = 3, q = 2 or p = 5, q = 3. Both cases can be treated in complete 
analogy to Example 2 above, and result in novel sum rules and combinatorial identities. 
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Recently we proved |T(J that exact trace formulae exist for one- dimensional square wells 
with an arbitrary number of potential steps inside. Following the methods outlined above, 
our results can be generalized immediately to obtain novel combinatorial identities for neck- 
laces with more than two types of beads. 
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j Wj Wj Vj n[Wj) 

1 CC C 2 1 

2 cn cn i 2 

3 nn n 2 i 



«K) /?K) t(wj) xK-) 

10 2 1 

2 3 2 

10 12 



Table I: List of the three cyclically non-equivalent binary necklaces of length 2 together 
with their primitives (w) and repetition indices (u). Some properties of the primitives, such 
as their lengths (n), number of TZ or C pairs (a), number of transmissions (/3), their weighted 
lengths (7) and their phases (x) are also listed. 



J 


Wj 


Wj 


lA i 


n(wj) 


a(wj) 


/3K-) 


7K-) 


x(wj 


1 


CCCC 


£ 


4 


1 


1 





2 


1 


2 


cccn 


cccn 


1 


4 


2 


2 


7 


4 


3 


coin 


cenn 


1 


4 


2 


2 


6 


5 


4 


cnai 


cn 


2 


2 





2 


3 


2 


5 


cjinn 


ennn 


1 


4 


2 


2 


5 


6 


6 


nnnn 


n 


4 


1 


1 





1 


2 



Table II: List of the six cyclically non-equivalent binary necklaces of length 4. The meaning 
of the columns is the same as in Table I. 
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